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OO \ We study the thermodynamics of massless (/)^-theory using screened perturbation theory. 

O 

. In this method, the perturbative expansion is reorganized by adding and subtracting a 



thermal mass term in the Lagrangian. We calculate the free energy through four loops 
expanding in a double power expansion m m/T and , where m is the thermal mass and 
, g is the coupHng constant. The expansion is truncated at order g'^ and the loop expansion 

. is shown to have better convergence properties than the weak-coupling expansion. The free 

■ ■ energy at order g^ involves the four-loop triangle sum-integral evaluated by Gynther, Laine, 

Qh| Schroder, Torrero, and Vuorinen using the methods developed by Arnold and Zhai. The 

Oh| evaluation of the free energy at order g"^ requires the evaluation of a nontrivial three-loop 

. sum-integral, which we calculate by the same methods. 

■ PACS numbers: ll.lO.Wx, 11.25.Db, ll.80.Fv, 12.38.Cy 

> 

OO 

^ : I. INTRODUCTION 

§ : n 

In recent years there has been significant progress in the understanding of thermal field theories in equilibrium MJ, 
4 13, U|- For example, the thermodynamic functions can be calculated as power series in the coupling constant g 
= at weak coupling and advanced calculational techniques have been developed in order to go beyond the first few 



' corrections. The pressure has been calculated through order g^ for massless 0'*-theory [S], massless QED 
and massless nonabehan gauge theories l3, [ill, Q ■ Very recently, the calculation frontier has been pushed to order g^ 
in massless ^''-theory by Gynther, Laine, Schroder, Torrero, and Vuorinen The calculation in Ref. Il3| involves 
the computation of complicated four-loop vacuum diagrams and was motivated by the corresponding problem in 
nonabehan gauge theories: There are three momentum scales - hard momenta of order T, soft momenta of order gT, 
and supersoft momenta of order g^T, which give contributions to the free energy. The contribution from the hard 
scale T to the free energy can be calculated as a power series in g^ using naive perturbation theory without resummed 
propagators. The order g^ is the first order at which all three momentum scales in QCD contribute to the free energy 
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and so it is important to calculate the full g^-term. Such a calculation involves the evaluation of four-loop vacuum 
diagrams in four dimensions. 

However, it is well known that the weak-coupling expansion is very sensitive to the renormalization scale, and it is 
furthermore convergent only if the coupling constant is tiny. The physical origin of this instability does not seem to 
be related to the magnetic mass problem in QCD, as it appears in 0'^-theory and QED as well. Rather, it seems to be 
associated with screening effects and quasiparticles. 

In recent years there have been large efforts to reorganize the perturbative series such that it has improved con- 
vergence properties. Several of these methods are variational in nature, in which the thermodynamic potential fl 
depends on one or more variational parameters rrii . The pressure and other thermodynamic quantities are then found 
by evaluating ft and its derivatives at the variational point where Sfl/Snii = 0. 

One of these methods is screened perturbation theory (SPT) which in the context of hot 0^-theory was introduced 
by Karsch, Patkos and Petreczky see also Refs. |l5l. Ila. Il7l|. In this approach, one introduces a single variational 
parameter which is added to and subtracted from the original Lagrangian. The added piece is kept as a part of the 
free Lagrangian and the subtracted piece is treated as an interaction. The parameter rn^ has a simple interpretation 
of a thermal mass and satisfies a variational equation. SPT has been applied to calculate the pressure to three- 



loop order 



18| and and the convergence properties of the successive approximations are dramatically improved as 



compared to the weak-coupling expansion. The mass parameter is of order g and so it might be reasonable to carry 
out an additional expansion of the Feynman diagrams in powers of m/T, and truncate at the appropriate order. This 
was done in Ref. 19(| and it was demonstrated that the double expansion in m/T and g converges quickly to the 
numerically exact result even for large values of the coupling. 

The generalization of SPT to gauge theories cannot simply be made by adding and subtracting a local mass term 
as this would violate gauge invariance. Instead one adds and subtracts to the Lagrangian a hard thermal loop (HTL) 
improvement term [20|. The free piece of the Lagrangian includes the HTL self-energies, while the remaining terms 
are treated as perturbations. Hard-thermal loop perturbation theory is a manifestly gauge invariant approach that 
can be used to calculate static as well as dynamic quantities in a systematic expansion. HTL perturbation theory has 

in an m/T expansion and the convergence 



been applied to calculate the pressure to two-loop order 



21 



properties of the successive approximations are again improved as compared to the weak-coupling expansion. 

Another variational method in which the propagator is a variational function was constructed by Luttinger and 
Ward 2^ and by Baym 2^1 for nonrelativistic fermions in the early 1960s. Later, it was generalized to relativistic 



quantum field theories by Cornwall, Jackiw and Tomboulis [23|. The approach is based on the fact that the thermo- 
dynamic potential SI can be written in terms of the two-particle irreducible (2PI) vacuum diagrams. The propagator 
D satisfies the variational equation Sfl/SD = 0. The 2PI effective action formalism is also referred to as ^-derivable 
approximations. 

Since the 2PI effective action formalism involves an effective propagator, a truncated calculation in the loop expansion 
or 1 /A^-expansion involves a selective resummation of diagrams from all orders of perturbation theory. This fact makes 
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renormalization of <i>-derivable approximations highly nontrivial. In recent years, there have been large efforts to prove 
renormalizability in the loop expansion, 1/A^-expansion, or the Hartree approximation, and in particular to prove that 
the counterterms are medium independent, i.e. independent of temperature and chemical potential 



29|, 3^,0,01. 



The second issue is that of gauge-fixing dependence. While the exact 2PI effective action is gauge independent at 
the stationary point, this property is often lost in approximations. The problem has been examined by Arrizabalaga 

3- In Ref. 

approximation. 



i nary 

m 



and Smit [33j as well as Carrington et al 



which is defined by the truncation of the action functional after n loops, has a gauge dependence that shows up 
at order g^". Furthermore, if the nth order solution to the gap equation is used to evaluate the complete effective 
action, the gauge dependence first shows up at order g'*". Explicit examples of the gauge dependence of the three-loop 
<i>-derivable approximation can be found in Ref. {s^. 

, 3^ and by Peshier Q| to calculate 
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The ^-derivable approach has been used by Blaizot, lancu, and Rebhan 
the thermodynamic quantities at the two- loop level in scalar field theory as well as in gauge theories. The calculations 
are based on the fact that the solution to the gap equation for the propagator for soft momenta is given by the HTL 



.e.t-e„e,.gie. Th»-.„„p calculations have been pe.fa.ed in scalar field theory hy Braa^ and Pe.itgirard Wi. and 
in QED in Ref. [33| using an m/T expansion similar to that employed in SET in Ref. The convergence of the 

successive approximations to the pressure is improved significantly compared to the weak-coupling expansion and 
the sensitivity to the renormalization scale is also reduced. In Ref. [4l|, the authors carried out a numerically exact 



three-loop calculation of the pressure in (^''-theory. Similarly, numerically exact two-loop calculations of the pressure 
in QED including an analysis of the gauge dependence of the results can be found in Ref. j^l • In these calculations 
no attempts to compare with the m/T expansions of Refs. were made. 

Finally, we mention other related resummation methods that have been applied in recent years, namely, the 2- 



particle point irreducable (2PPI) method 
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These methods 



4J| as well as the linear delta-expansion 
are also variational in spirit. Moreover, it has been shown that they correctly predict a second-order phase transition 
when applied to 0*-theory. In the case of the linear delta-expansion, the successive approximations of e.g. the pressure 
are remarkably stable as compared to the weak-coupling expansion. 

The article is organized as follows. In Sec. II, we briefiy discuss the systematics of screened perturbation theory. 
In Sec. Ill, we calculate the pressure to four-loop order in a double expansion in m/T and . In Sec. IV, we discuss 
different gap equations that are used to determine the mass parameter in screened perturbation theory. We also 
present our numerical results and compare them with the weak-coupling expansion. In Sec. V, we summarize. In 
Appendix A and B, we list the sum-integrals and the integrals that we need. In Appendix C, we discuss the m/T 
expansion of typical sum-integrals that appear in the calculation. In Appendix D, we calculate explicitly a new 
three-loop sum-integral that contribute to order g'^ in the m/T expansion. 
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II. SCREENED PERTURBATION THEORY 



The Lagrangian density for a massless scalar field with a (p* interaction is 

£ = i5^09^0-|^04 + A£, (1) 

where g is the coupHng constant and A£ includes counterterms. Renormalizability guarantees that AC is of the form 

1 . , , 1 



AC = ^AZ d^4>d''(t> - f^^9^<t>^ • (2) 

Screened perturbation theory, which was introduced in thermal field theory by Karsch, Patkos and Petreczky [l^ . 
is simply a reorganization of the perturbation series for thermal field theory. It can^be made more systematic by using 
a framework called "optimized perturbation theory" that Chiku and Hatsuda 
broken scalar field theory. The Lagrangian density is written as 



401 have applied to a spontaneously 



CsPT = -fo + - o ("^ - O 4^ " ttA + A£ + A£sPT . (3) 



Here, £o is the vacuum energy density term, and we have added and subtracted mass terms. If we set £q — Q and 
m\ ~ m? , we recover the original Lagrangian Eq. ([1]). Screened perturbation theory is defined by taking to be 
of order unity and m\ to of order g^, expanding systematically in powers of and setting = m? at the end of 
the calculation. This defines a reorganization of the perturbative series in which the expansion is about the free field 
theory defined by 



The interacting term is 



£int = \ml(tP - ^<P^ + A£ + A£sPT . (5) 



Screened perturbation theory generates new ultraviolet divergences, but they can be cancelled by the additional 
counterterm in A£spt- If we use dimensional regularization and minimal subtraction, the coefficients of these operators 
are polynomials in g"^ and {w? — ml). The counterterm A£ is 

A/: = -^c^'. (6) 

The additional counterterms required to remove the new divergences are 

A/:sPT = -A^o - I {Am^ ~ Ami) ■ (7) 

Several terms in the power series expansions of the counterterms are known from previous calculations at zero temper- 
ature. The counterterms Ag^ and Am^ are known to order a^, where a — g^/(47r)^ j^l- We will need the coupling 
constant counterterm to next-to- leading order in a: 



Ag' 



3 / 9 17 , 



2 



9' ■ (8) 
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We need the mass counterterms Am^ and Am^ to next-to-leading in a: 



Ato^ 
Atoi 



2e 
1 



1 5 
1 5 



(9) 
(10) 



2e ' V2e2 2Ae, 

The counterterm for A5o has been calculated to order 5l||. We will need its expansion only to first order in a and 
second order in m^: 



(47r)2Afo 



" 1 


1 




■ 1 




2 2 1 


■ 1 
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8e2" 
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III. FREE ENERGY TO FOUR LOOPS 



In this section, we calculate the m/T expansions of the pressure to four loops in screened perturbation theory. In 
performing the truncation, m is treated as a quantity that is 0{g) and we include all terms which contribute to order 



A. One-loop free energy 

The free energy at leading order in is 

^0 = + ^Oa + Aofo , (12) 
where Aq^^o is the term of order in the vacuum energy counterterm Eq. (fTTj) . 



o 



Figure 1: One-loop vacuum diagram. 
The expression for diagram J^oa in Fig. [T] is 

= log [P^ + m^] , (13) 



where the symbol ^ is defined in Appendix [Al 

Treating m as 0{gT) and including all terms which contribute through 0{g'^), we obtain 



■^oa - + im^Xi + ir/^ - + im% , (14) 
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where the sum-integrals 2q and X„ are defined in Appendix [A] and the integral Iq is defined in Appendix [BJ In 
Appendix C, we illustrate the m/T expansion of simple one-loop sum-integrals such as the one appearing in Eq. lfT3|) . 
We also note that most of the multiloop diagrams are products of simple one-loop sum-integrals. 

The term X2 is logarithmically divergent and the pole in e is cancelled by the zeroth-order term AqSq in Eq. ifTTj) . 
The final result for the truncated one-loop free energy is 



^0 



where m = and L — log . 

2.7:1 47ri 



90 



15 



1 - 15m^ + 60to^ + 45(L -t- -fE)m'' - yC(3)TO'' 



(15) 



B. Two-loop free energy 



The contribution to the free energy at two loops is given by 



Ti = .Fia + .Fib + Ai£o + AiTO^ ^ 



(16) 



where Ai^o and Aim'^ are the vacuum and mass counterterms of order 5^, respectively. The expressions for the 



00 o 



la 



Figure 2: Two-loop vacuum diagrams. The cross denotes a mass insertion. 



diagrams jFia and Tih in Fig. [2] are 



1 



1 



Fib = ^2^1 

Expanding the sum-integrals in Eqs. (fT7|) and lfT8|) to order 0{g'^) yields 
1 



Fla 



ll + 2TI{Ii - 2m^l{l2 + T^ii - 2m^IiTl2 + 2to''XiX3 + m^I^ + 2m^TI{Ii 



1 



■Fib = ^-jj^i 



(17) 
(18) 

(19) 
(20) 



where the integral /„ is defined in Appendix [Bl 
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The poles in e in Eqs. (fT9| and l(20|) are cancelled by the counterterms in Eq. (fT6|) . The final result for the two-loop 
contribution to the free energy is 



c / 1 \ 45 

- - 15m - 15 (L + 7B - 3) + 90(L + ^ieW + 45 [ (i + 7is)' + —((3) J - y C(3)m^ 



90 



15m? 



3 

1 - 6m - 6(L + 7ij)m2 ^ 7:C,{'i)m^ 



(21) 



Note that we here and in the following have pulled out a factor of ^idcai = — tt^T'^/QO 



C. Three-loop free energy 



The contribution to the free energy at three loops is 



^2 — + ^2h + ^2c + ^2d + ^2^0 

^:^Ai5^ + 5^Aim?, 



A^m H ^/A.™^\ II I 1 A ™^ 

2 (am2)' 



dm? 



(Aim^ 



\ 9m2 Sm^ 



Aim^ 



(22) 



5 "'1 

where we have included all necessary counterterms. The expressions for the diagrams JF2a, ^2b, .?^2c, and !F2d in Fig. [3] 



ooo 




T2. 



^2h 



00 o 



^2c 



Figure 3: Three-loop vacuum diagrams. 



•?^2a 
•^2b 
^2c 
•^2d 



m 
1 



2^2 



_1 4y/ 1 1 

48^ l^PQR P2 + m2 Q2 _^ m2 i?2 _^ m2 (P + Q + i?)2 + 



2\2 " 



(23) 
(24) 
(25) 
(26) 
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Expanding in powers of to the appropriate order gives ^ 
1 



^2, 



16^ 



2b 



2m2r%/i/2 - 2m^Iill - Am^Thlils - 2m^TIiI^ + m^Thll + 2m^Tl2liTz , (27) 
1 . _ . „ , ^28) 



48^ 



T 1 2 2 

^2c = -^g 



2d 



Tl2+l2-2m^l3 



, (29) 
(30) 



where Isun, and Ibaii, and I^aii ^^'s defined in Appendix [AJ and /bail is defined in Appendix [BJ 
The poles in e in Eqs. (|27|) - (|30)) are cancelled by the counterterms in Eq. (|22l) . 
The final result for the three-loop contribution to the free energy is 



^2 = -- 



t:^T^ 5 1 



90 8 TO 



42to2 5 + 7L + 37_e - Slogm - 81og2 - 4 



C(-i) 



268(L + 7b) - 48(L + 7^)^ + 



(34 + 127b) 



+12 



C"(~l) . „. 9tt^ 



+ 7B (17 - 217b) + 34 + — - 4871 - C(3) - eC^.n to' 



C(-l) ■ ■ ' 2 

+ ( 89 + 120(L + 7b) + [18(L + 7^)]^ + 15C(3) ) to"* 



TT^T* 15 mi 

H 

90 2 TO 

TT^T^ 45 mf 
90 2"^ 



+ 2 (L + 7B - 3) TO - 18(L + 7B)m2 _ (^-^^ (L + -fsf + C(3)) + y C(3)to4 
1 + 2(L + 7b)to-C(3)to^ 



(31) 



Here C^^ii ~ 48.7976 is the numerical constant in I^aii 



^ Notice that the term T/iTsun in T2b in Eq. eqref2b vanishes. However, we include this term because it gives rise to a finite 
term at four loops when renormalizing the coupling constant g. 
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D. Four- loop free energy 



The contributions to the free energy at four loops are 



(pm'^y \ om'^ om'^ / g'^ \ 9 9 9 

+ 1 (T^k + T^k) (Ai-^)^ + ^A^n.? + (A,..^) (A.ml) + 1^ (A,,^) (A.m^) 



2 \(9m^)^ [dm'^Y J m\ m\dm? ^ ^ 9^ dm? 

+ 2-^ Aim? + — ^ + — ^ + — ^ + — ^ Aim2 . 32 
mf J \ om'^ om^ om^ om-' J 

Note that some of the terms first contribute at order or higher. For example, the vacuum counterterm A^Sq first 
contributes at order m'^a'^ ~ g^. 

The expressions for the diagrams J-sa-J^sj , in Fig. |4] are 

2 / „ -, \ 2 



^3a 32-9' (^p2 +^2) (^(Q2+^2)2) 



3e 



^2-)3 ' 



(33) 



24'^ ypQR{P^+m^)^Q'^+m^R^+m^{P + Q + R)^ + m^%S'^+m^ ^ ' 

^ = 1 6v/ ^ 1 I I I 1 



(38) 



•^3g - ~Y^g '"1^^^^ (p2 + ^2)2 Q2 + ^2 ^2 + ^2 (p + g + p)2 + ^2 (^9) 
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oooo 





© ooo 



OOO 




^3f 



•^3g 



oo oo 



^3h 



o 



Figure 4: Four-loop vacuum diagrams. 
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Expanding the sum-integrals in powers of to the appropriate order gives 



■^3a — 

^3c = 

^3d = 

^Za — 

Tzi = 

^3g = 

^3h — 

•^3j - 



32 

-2m^T^Iilil2 + 472/1/2X1X2 - 'im^Thlill + 2X^/2/2X2 - 4m2T/2Xi2X3] , 

48 

+3TIiI^l3 - ?,m^T^llhl2 + im^Thlll^ + Sm^'X/aXiXl] , 

.96 [(Xi + T/i - 777,2X2) r3/(,,ii + (Xi + T/i)X(,,ii + 3X2/1/2X1X2 - m^Tl2l2lsnn 



24 



Q2 + (2/d)q2 



+ 3X3/2/2X2 + 2X1 (T/i - 7„2^,^^ ^^^^^^ ^ 

^.9' [n(P)]' + X^/tHanglc + 6r/l^ [n(P)]' + 3T3X2/ban 



1 



.9^7772 [x2/22Xi + T3/i/22 + 2X/2X1X2 - 277727/2X1 + X1X22 + 2x2/1/2X2 - m^T^lll2 
+TIili -Am^T I2I1I3] , 

-^g^ml [7/3X2 + 2x2/3/1X1 + 73/3/2 - 277727/3X1X2 +X3X12 - 2727772/3/1X2 + 27/1X1X3 
+2777^^7/3X1X3 + 777^7/3X22] , 

[73/;,„ + XLn + +372/1/2X2 + 2 (7/1 - 77727/2) + 

Jo. 



12 



4 2 
.9 ™i 



-52^4 [7/3x1 + 72/1/3 + 7/1x3 - 77^27x2/3 +X1X3 + 777^7x3/3] 

^52^4 [j,2^2 ^ 27/2X2 +X22 - 477727/2X3] , 

-i7n?[7/3+X3] , 
D 



(43) 
(44) 

(45) 
(46) 

(47) 

(48) 
(49) 

(50) 
(51) 
(52) 



where the self-energy 11(7) is defined in Eq. l|D.2p and the integrals /^^^jj and /triangle are defined in Appendix B. The 
poles in Eqs. l(43l) - l(52|) are cancelled by the counterterms in Eq. l(32l) . The final result for the four-loop contribution 
to the free energy is 
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90 288^ 



1 + 18( llL + 37E-6-161og2-161ogm-8^^^ Im^ 



/ Q1Q8 6456 

- ( 1236 + 108q^^„gi, + 36C;^ii + 28871 - -^1e + 4507| - -^L + 4327bL + 648^^ + 1357r2 



547r^Ci,_ie - 2167r^7B + (2100 - 727^ + 1728X) 



C(-i) 



432 



C(-i) 



432 (7B + 2L) 



C(-3) 
C(-3) 



+360 



C"(-i) 
C(-i) 



+ 1728 log 2 + 21677^ log2 + 432(4 - tt^) logm - 4534C(3) m 



+ ^ ^3742 - 288(7/ - 48C(,aii - 8O6471 - 60727b - 25447| - 3904L - 18727eL 



-2184L2 + qoOtt^ + (1808 + 18247/=; + 2496L) 



C'(-i) 



- 288 



C"(-i) 



+ 26887b log 2 



C(-i) C(-i) 

+4992ilog2 + 4992(7/5 + L) logrfi - 23047/jlog7r + 23041og^(27r) - 15C(3) \ 
84L + 36jE- 96 log m - 36 - 96 log 2 - 48 \ ' _ / | 



90 16 m3 ' ' '"C(-l) 

+2 f 48(L + 7/;)2 - 268(L + 7/5) - 7/=;(17 - 217/5) + 4871 - 34 - — 

- fpiy (34 + 127/5) - 12^p^ + C(3) + 6C;,u) - 3 (^89 + 120(i + 7i=;) + [18(L + 7i=;)]' 



+ 



+15C(3) jm^ 
TT^T" 15 ami 



90 8 m3 
_7r2r4 15 mf r 
90~Trf!3 

where the constants are 



45 



1 + 18(L + 7/5)m2 ^ [24(i + 7^)^ + 2C(3)]m3 - y C(3)m^ 



1 + 2C(3)7fi^ 



(53) 



CU = 48.7976, 

^"riangle = —25.7055 , 

<^triangle = 28.9250, 

Ci = -38.5309. 



(54) 
(55) 
(56) 
(57) 



There are a couple of calculational details that are worth while pointing out. The contribution arising from diagram 
J-'zd when all momenta are hard reads 



^(hhhh) _ 



= ^[n(P)]^ 



(58) 



This term can be combined with the term arising from the counterterm J^2h^ig^ = —9^^ha.n^i9^/^8 and gives 

3 



i 



(59) 
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This particular combination was first calculated by Gynther et al. [12] using the methods of Arnold and Zhai. Similarly, 
we combine the term from J-^d with the term T/iXgun from J^2bAi5^, which gives 

gi^{p(^)]^-(4^P(P)]}. (60) 

We calculate this sum-integral in Appendix D. Finally, the term from JF2bAim^ which involves Isun can be combined 
with the term — m^/22'22'sun arising from to give 

Since Xsun vanishes at order e° and the term inside the paranthesis is finite, the particular combination (|6T|l vanishes 
in the limit e ^ 0. 



E. Pressure to four loops 



The pressure V is given by — JF. The contributions to the pressure of zeroth, first, second order, and third order 
in are given by Eqs. (flSl) . i(2T|) . (|3T|) . and (|53l) . respectively. Adding these contributions and setting Sq = and 
nil — "^^) obtain approximations to the pressure in screened perturbation theory which are accurate to O (g^). 
The one-loop approximation to the pressure is 



Vo = T^idoal 



1 - 15m2 + 60to3 + Abm^{L + 7_e) - ^((3) 



(62) 



where T^idcai — 7r^T'*/90 is the pressure of an ideal gas of massless particles. 

The two-loop approximation to the pressure is obtained by adding Eq. l(2T|) with nil = w?: 



Vo+i 



'Pidcais 1 — -^oi + 15rha + 15fh {L + je — 3)a — 30m 



l + 3(L + 7£;)a 



(L + 7B)+((i + 7B)' + ^C(3))a 



45 

-C(3)m'^a + 15C(3)7 



(63) 



The three- loop approximation to the pressure is obtained by adding Eq. ((311) with m\ =w?: 

Vq+1+2 



5 2 5 
■Pidcal< 1 + -TT^a ~ -a 
8m 4 



59 15^ 5 ,C'(-1) 5C'(-3)\ 2 

i2+T^+4^-+^cRy"2^r 



15 

7 

2 

Ci-i) 
C(-i) 

15 



1- ( 5-|-37j;-|-7L-81ogm-81og2-4 ^,f ^} \a 



C(-i) 



a + -m^ ( 268(L + 7e) - 48(L -f 7^) 



(34+ 127£;) + 12 



r(-i) 

C(-i) 



7£;(17 - 217b) + 34 + — - 4871 - C(3) - 6C; 



ball 



1 - 6(L + 7£)a - ( 89 -f 120(L + je) + [18(L + 7^)]^ + 15C(3) ) a' 



+45m^ ( (L + ^E? + ^C(3)) a - ^m^aC(3) - ^rn'm 



(64) 
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The four- loop approximation to the pressure is obtained by adding Eq. ((53l) to Eq. with ml — vn?: 

-3 15 1 



^0+1+2+3 _ 5 a"- 

■Pidoai ^ 288 m3 16 TO 



+ 1 (^16 log m + 6 - 37is - IIL + S^p^ + 16 log 2 ) 



« _ ( 3£ _ + - 2^^) + ^ ( 1236 + 36CL„ + 108^'^, 



+28871 



9198 



-7b + 45071; 



6456 



L + 4327BL + 648L^ + 1357r^ - ^^t^"^ C^-.^^^y^ 



-21&TT'^"fE + (2100 - 727B + 1728L) 



C'(-i) 
C(-i) 



432 



C(-i) 



432(7B + 2L 



C(-3) 
C(-3) 



45 



+360 + 1728 log 2 + 2167r2 log 2 + 432(4 - tt^) log m - 4534C(3) ) 
2 



a-- yli + i^E + lL-A 



^ -81og2- 81ogmj ~ ^( 3742 - 288C/ 



-48C{,^n - 806471 - 60727^ - 25447| - 3904L - 1872-fEL ~ 21ML^ + 9007r2 

C'{-1) C"{-1) 
-(1808 + 18247b + 2496L)-^^^ - 288-^^^-^ + 26887B log 2 + 4992L log 2 

-4992 (7B + L) log TO - 23047B log 7r + 2304 log2(27r) - 15C(3) 
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4 
135 



TO^ 



1 - 3(L + 7B)a + ^ ( 89 + 120(L + 7^) + [18(X + 7^)]" + 15C(3) ) a' 



C(3)m5a 



(65) 



The final result for the pressure is given by Eq. (|65|) . If we use the weak-coupling expansion for the mass parameter, 
TO^ = a/6, our result reduces to the weak-coupling expansion result through order . ^ Inserting to^ into Eq. l(65|) . 
we obtain 



5\/6 



V = "Pidoai"! 1 ~ 4" ■ 3 



,3/2 



15 



log 



AttT 



C4 



45 
'T 



log 



2 fJ- 



1 /269 



15\/6 



2. oC^(-l) , 
AttT 3 V 45 C(-l) C(-3) 



log 



1 



(66) 



where the constants C4-C6 are 

59 1 4C'(-1) 2C'(-3) 



C5 



45 3^^ 3"C(-1) 3 C(-3) ' 

5 1 2 2 2 C'(-1) 
77 + 7;lE - - log ; 



3 C(-l) 



103 1 



511 



-(4 - TT^) log - + — + — Cball - cQrianglo " T^C'trianglc + Ti^l ~ 



25 



54 18 



6 



12 



180'^^ 36^^ 



y, „ .175 1 

+7r'^ log2 + 7b 

6 ^ I 54 9 '-^ 



C'(-l) 2 /C'(-1)A 5C"(-1) 2 C'(-3) 2267 



C(-i) 3VC(-i) 



+ - 



9 C(-l) 3^^"C(-3) 



324 



24 



C(3). 



(67) 
(68) 

(69) 



It is important to point out that we have only calculated part of the y'^-term in the weak-coupling expansion. See the 
discussion in Sec. [V] 
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The numerical values of C^-Ce are 

C4 = 1.09775, (70) 

C5 = -0.0273205, (71) 

Ce = -6.59363. (72) 

Gynther et al. 



13] have calculated the pressure for an 0(A^)-symmetric theory at weak coupling through order 

using effective field theory methods. Our result agrees with theirs for = 1. 

Using the renormalization group equation for the running coupling constant to next-to-leading order, 

da 9 17 , , , 

= 3a^ - -a3 , (73) 

it is straightforward to verify that the result ((66| is independent of the renormalization scale /k through order log (7. 



IV. GAP EQUATIONS AND NUMERICAL RESULTS 

The mass parameter m in screened perturbation theory is completely arbitrary. In order to complete a calculation 
using SPT, we need a prescription for the mass parameter m as a function of g and T. One of the complications 
which arises from the ultraviolet divergences is that the parameters £0, , m\, and g^ are all running parameters 
that depend on the renormalization scale /x. 

The prescription of Karsch, Patkos, and Petreczky for m,(T) is the solution to the one-loop gap equation: 



Ji{f3m^)T^~ (21og-^ + l ) 



(74) 



where /i* is the renormalization scale and Ji(/3m) is the function 

Their choice for the scale was /i* = T. In the weak-coupling limit, the solution to l(74|) is m, = g(/i*)T/V24. The gap 
equation l(74|) is the renormalized version of the following equation 



There are many possibilities for generalizing ([74| to higher orders in g. We will consider three different possibilities 
in the following. 

1. Debye mass 

One class of possibilities is to identify to* with some physical mass in the system. The simplest choice is the Debye 
mass Tod defined by the location of the pole in the static propagator: 

p2 + + I](0,j3) = 0, p^^-ml. (77) 
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The Debye mass is a well defined quantity in scalar field theory and abelian gauge theories at any order in perturbation 
theory. However, in nonabelian gauge theories, it is plagued by infrared divergences beyond leading order js^ . 



2. Tadpole mass 

The tadpole mass is another generalization of Eq. l(74|) to higher loops. It can be calculated by taking the partial 
derivative of the free energy T with respect to before setting rrii ~ m: 



2 2 
m-t = g 



(78) 



From this equation, we see that TOj is proportional to the expectation value {4>^). The tadpole mass is well defined at 
all orders in scalar field theory, but the generaHzation to gauge theories is problematic. The natural replacement of 



a2\ 



would be {Af^A^j), which is a gauge- variant quantity. 



3. Variational mass 

There is another class of prescriptions that is variational in spirit. The results of SPT would be independent of m 
if they were calculated to all orders. This suggests choosing m to minimize the dependence of some physical quantity 
on m. The variational mass is defined by minimizing the free energy; 

The variational mass has the benefit that it is well defined at all orders in perturbation theory and can easily be 
generalized to gauge theories. 



4- Comparison 

At one loop, the three different prescriptions give the same gap equation, Eq. ((74|) . Moreover, it turns out that 



the two-loop tadpole mass coincides with the one- loop tadpole mass However, at two loops the screening and 

variational masses are ill-behaved [l^. The screening mass solution ceases to exist beyond g ^ 2.6 and the variational 
gap equation only has solutions in the vicinity oi g — for some values of L. In the following, we therefore restrict 
ourselves to the tadpole gap equation. 
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5. Tadpole gap equation through three loops 



At one loop, the renormalized gap equation follows from Eq. (fT5| upon differentiation with respect to w? and can 
be written as 







a 



1 - 6m - 6m2 (L + 7_e) + -Q{i)rh'^ 



(80) 



At two loops, the renormalized gap equation follows from differentiating the sum of Eqs. (|T5l) and (|2T|) with respect 
to m, and setting mi = m. It can be written in the form 



= m^ 



12m 6 



-22 

-m a 



1 + a(3 - 7B - L) 
5 



1 ^ 

—ma 

2 



1 - 3a(7£; + L) 



, C(3) 
(7b + ^) + 



1 



(81) 



At three loops, the renormalized gap equation follows from differentiating the sum of Eqs. (flSl) . l(2T|) . and (|3T|) and 
setting mi = m. This yields 







1 a' 



1 + a 



7^ 4C'(-1) 8, „ 8, / 
l-7.-3^+3^ + 3log2 + -logm 



l-a(L + 7B-3) + a2 



2(L + -lEf - + 271 - ^(i + 7b) - ^7^(17 - 217^) 



3^ 
16 



i7C'(-i) 1 C'(-i) iC"(-i) 



1 



1 



12 C(-l) 2^"" C(-l) 2 C(-l) 



^ + ^C(3) + 



+ -?7ia< 1 — 2a{L + 7^) + a 



9(i + 7i=;)' + y(^ + 7i^) + § + ^C(3) 



^m3a\(3) 



(82) 



6. Numerical results 



The two-loop SPT-improved approximation to the pressure is obtained by inserting the solution to the one-loop 
gap equation l(80|) into the two-loop pressure (|63|) . In Fig. [SK we show the various truncations to the two- loop SPT 
improved approximation to the V /V ideal as a function of g{2TTT). We notice that the various truncations converge 
quickly. The order-g^ to order-g^ results are almost indistinguishable and essentially equal to the exact numerical 
two-loop result in Ref. In the three- loop case, we insert the solution to the two- loop gap equation (HJ) into the 
three-loop pressure l(64|) . In Fig.jSb, we show the various truncations to the three-loop SPT-improved approximation 
to V /V ideal as a function of g{2TrT) . The three-loop result also converges to the exact numerical three-loop result, 
albeit not as fast as in the two-loop case. At four loops, we insert the solution to the three-loop gap equation l(82|) into 
the four- loop pressure (|65| . In Fig. [5b, we show the various truncations to the four- loop SPT-improved approximation 
to 7^/7^idcai as a function of g{2'KT) . Although we cannot compare our successive approximations with a numerically 
exact 4-loop result for the pressure, we expect them to converge reasonably fast. Based on the experience with the two- 
and three-loop approximations, we expect that the g^-truncation provides a good approximation to the numerically 
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g{2TTT) 

Figure 5: (a) 2-loop pressure, (b) 3-loop pressure, (c) 4-loop pressure, (d) weak-coupling expansion of the pressure, 
all normalized to T^idoai. 

exact result. Clearly, however, only a calculation through can settle this issue firmly. In Fig. [Hi, we show the 
weak-coupling expansion of 7^/7^idcai to orders g"^, g^, g^, g^, and g^ as a function of g{2'KT) for comparison. Note 
that the results to order are identical in SPT and in the weak-coupling expansion since there is no m-dependence 
at this order. 

In FiglHK, we show the two, three- and four- loop pressure through order g"^ normaHzed to V /V ideal as a function of 
g{2nT) . In Fig. [Bb, we show the weak-coupHng expansion of "P/Pidcai to orders g^, g^, g'^, g^, and g^ as a function of 
g{2TrT) for comparison. The successive approximations using screened perturbation theory have better convergence 
properties than the weak-coupHng results. The improved stability is partly due to the fact that we are using a thermal 
mass determined by a gap equation and not by the perturbative value for the Debye mass. 
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Figure 6: (a) Pressure normalized to 'Pideai through for various loop orders, (b) weak-coupHng pressure at various 
orders of g. 

V. SUMMARY AND OUTLOOK 



In this paper, we have calculated the pressure of massless scalar field theory to four loops using screened perturbation 
theory expanding in a double expansion in powers of g^ and m/T. Treating m as 0{gT), we truncated our expansion 
at order 5^. The expansion required the evaluation of a new nontrivial three- loop diagram, where we evaluated the 
sum-integral using the techniques developed in Ref. We have seen that the successive approximations are more 
stable than the weak-coupling expansion. In particular, it is interesting to note that the four-loop curve lies between 
the two-loop curve and the three-loop curve. The apparent improved convergence seemed to be linked to the fact that 
SPT basically is an expansion about an ideal gas of massive particles instead of an expansion about an ideal gas of 
massless particles which is the case for the weak-coupling expansion. 

Using the weak-coupling value for the mass parameter m, our result reduces to the weak-coupling result for the 
pressure through g^. In particular, we have reproduced the pressure at weak coupling for = 1 obtained by Gynther 
et al [1^ • Using effective-field theory methods , the authors in Ref. l3| have calculated the hard and soft contributions 
to the pressure through order g^ separately. It appears that the convergence properties in the hard sector are better 
than in the soft sector even for moderate values of the coupling. 

We have mentioned that our result only includes part of the full ^''-term in the weak coupling expansion. This is 
straightforward to see if one uses the effective-field theory approach developed in j^. The contributions to the free 
energy comes from the two momentum scales T and gT. The contribution from the hard scale T can be calculated 
by evaluating the sum-integrals with bare propagators and so is therefore a series in g^ starting at order g^. The 
contribution to the free energy from the soft scale gT can be calculated using an effective Euclidean three-dimensional 
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field theory whose coefficients depend on g and T . This contribution to the free energy is a series in g starting at g'^ . 
The contributions to the free energy that are odd in powers in g is therefore entirely coming from three-dimensional 
vacuum diagrams and power-counting tells you immediately that part of the 5^ term is arising from the five-loop 
vacuum diagrams. Our four- loop calculation therefore agrees with the weak-coupling expansion through order g^ . 

In order to evaluate the free energy to order g^, we must determine all the coefficients in the effective theory to 
sufficiently high order in g. The only nontrivial calculation that is required is to determine the mass parameter in the 
effective theory to order 5®. This involves the expression for the diagram calculated in Appendix D i.e, the sum- integral 



/ = 



5f^{[n(P)P-(4^n(P)}, (83) 



The evaluation of the free energy to order g'^ is in progress 



3- 
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Appendix A: SUM-INTEGRALS 

In the imaginary-time formalism for thermal field theory, the 4- momentum P = (Po,p) is Euclidean with = 
Pq + P^- The Euclidean energy po has discrete values: Po = 2n7rT for bosons, where n is an integer. Loop diagrams 
involve sums over Pq and integrals over p. With dimensional regularization, the integral is generalized to d = 3 — 2e 
spatial dimensions. We define the dimensionally regularized sum- integral by 

where 3 — 2e is the dimension of space and [i is an arbitrary momentum scale. The factor [e^ I^-kY is introduced so 
that, after minimal subtraction of the poles in e due to ultraviolet divergences, [i coincides with the renormalization 
scale of the MS renormalization scheme. 



1. One-loop sum-integrals 

The massless one-loop sum-integral is given by 



'P 



....... C(2n - 3 + 2e) r(n - I + e) .^^j..4-2n-2. 2) 

8^2 r(i)r(n) ^^""^^ ' ^^-^^ 
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where C{x) is Riemann's zeta function. Specifically, we need the sum-integrals 



T' 



^3 



logP 

p 



45 



12 UttT/ 



C'(-l)\ / 



C(-i) 



1 



2e 



1 



(47r)2 V47rr 
(i^[2C(3) + 0(e)] 



27B 



- 471 e +0 (e^ 



(A.3) 
(A.4) 

(A.5) 
(A.6) 



2. Two-loop sum-integrals 



We need two two-loop sum-integral that are Hsted below: 

1 



4 



Is. 



p2 + (2/d)p2 



4 

JF 



0{e), 



pqP'Q^{P + Q? 4(4^) 



ttY \AttTJ 



1 



10 



^+hj+47£; - + :;77+o + -irlE + 47I; - 871 + 0{e) 



(A.7) 
(A.8) 



V6 ' J e 36 2 3 

The setting-sun sum- integral was first calculated by Arnold and Zhai in Ref. [lJ|, while Eq. I|A.8|) was calculated in 
Ref. 3- 



3. Three-loop sum-integrals 



We need the following three-loop sum- integrals: 
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|4 V47rTy 
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46 - 87B - I67I; - 10471 - 247£; log(27r) + 241og^(27r) 
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1671 



C'(-i) 



+ Ci + 0{e) 



(A.9) 



(A.IO) 



457r2 



(A.ll) 



where C^^n 



C(-i) C(-i) ' ^"'"C(-i) 

48.7976 and C/ = —38.5309. The massless basketball sum-integral was first calculated in Ref. [IQ], and 
X^jjji in Ref. |13|] . The expression for the sum- integral Eq. (|A.lip appears here for the first time and is calculated in in 
appendix D. 



4. Four-loop sum-integrals 

We also need a single four-loop sum- integral which was calculated in Ref. [l^: 



^^{[n(P)]' 



[Tl{P)f = 



16(47r) 



where Q,;^^ 1, 



27£ + 4^^^ ) (2L + 7b) + Qnanglc 



C(-3) 



-25.7055 and C'f,;,^!, 



28.9250. 
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3 ^ C(-i) 



4L ) - + (2L + ^e) 



512(47r)2 
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+ ^L + 47b + C( 



triangle 



+ 0{e), (A.12) 



Appendix B: THREE-DIMENSIONAL INTEGRALS 



Dimensional regularization can be used to regularize both the ultraviolet divergences and infrared divergences in 
3-dimensional integrals over momenta. The spatial dimension is generalized to d = 3 — 2e dimensions. Integrals are 
evaluated at a value of d for which they converge and then analytically continued to c? = 3. We use the integration 
measure 

(B.l) 



eV 
47r 



(2^) 



3-2e 
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1. One-loop integrals 



The one-loop integral is given by 



Specifically, we need 



1 (e7«^2).£|L_l+i)^3-2„-2. 



Stt^'^ ^ ' r(i)r(n 



/l = 
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log(p^ + m^) 

p 

677 V2m/ 

m / fi \ '^'^ 
Itt V2m/ 

87rm \2m/ 

^— f— V 

3277771-^ \2m/ 



l + 2e+(4+— le^ + O (e^^ 



1 + 2e + — + O (e3) 
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(B.4) 
(B.5) 
(B.6) 



2. Three-loop integrals 



We need two three- loop integrals: 

111 1 
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The massive basketball was calculated in Ref. to order e°, and to order e in Ref. [sJl- The other 3-loop integral is 
obtained by differentianting the massive basketball with respect to the mass m. 



3. Four-loop integrals 



We need a single four-loop integral, namely the triangle integral. This integral was calculated in Ref. j3| and reads 

111111 



^triangle — 
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Appendix C: m/T EXPANSIONS 



In this appendix, we list the m/T expansions of the sum-integrals we need. The sum-integrals include sums over 
the Matsubara frequencies Po = 27rnT and integrals over the three- momentum p. In the sum- integrals, two important 
mass scales appear. These are the hard scale 2ttT and the soft scale m. The soft scale m is of order gT and at weak 
coupling this scale is well separated from the hard scale, m ^ 2ttT . We can therefore expand the sum-integrals as a 
Taylor series in powers of m/T. 

First consider the simple one- loop sum-integral appearing in the expression for the one-loop free energy in Eq. ifTS]) : 



•^Oa = i^log[p2+^2] 



1 /-(h) 1 

= log + + 2% + ' (^-^^ 

where the superscripts (h) and (s) denote the hard and soft contributions, respectively. In the hard region, the 
momentum P is of order T and so we can expand in powers of m? /P"^ . This yields 

log + „,2] ^ ^ p2 ^ _1_ _ 1 „^4^ + • • • (C.2) 

The contribution from soft momenta is given by the po = mode alone and and reads 

^^\og [p2+m2] = T J^logip^ + m^) . (C.3) 

The other simple one-loop sum-integrals are expanded in a similar manner. 
We next consider the massive basketball diagram in Eq. l(24l) : 



Xban(m2) + ^2)(Q2 + ,^2)(^2 + ^2)[(p + Q + i^)2 + ^2] " (^-4) 

Eq. (|C.4p involves three sum- integrals and so receives contributions from four momentum regions: (hhh), (hhs), (hss), 
and (sss). In the first case, where all the loop momenta are hard, we can expand the sum-integral in powers of m^. 
This yields 



-(hhh) 
ball 



fpQR P'Q'RHP + Q + R? fpQR P^Q^R'iP + Q + Rf 
When two momenta are hard and one is soft, the contribution reads 

I^'^'^'Hm') ^ AT f - y/ - - - 

^ ' JpP^ + m^%i,Q^ + m^R^ + m^p + Q + R)^+m'' 



= AT 



When one momentum is hard and two are soft, the contribution is given by 



(C.6) 



^all (m j - y^^^2+^2^2 + ^2^^^2 + 



m^ (p + q + R)^ + 

/ ^2 :^2.2 :^2 a • (c.7) 
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Finally, when all momenta are soft, the contribution is given by the massive basketball diagram /bail in three dimensions: 

^baii ("^ ) ^ J^^^ p2 _|_ g,2 _|_ ^2 J.2 + 777,^ (p + q + r)^ + ■ i^-^) 

The basketball diagram with a single mass insertion, 2'^aii(™^)' '^^^ calculated by differentiating the massive 
basketball diagram with respect to m^. This yields 

^, ( 2. ^ 1 1 1 1 

balll"^ ) ^^^^ ^p2 ^ ^2)2 g2 + ^2 ii;2 + m2 (P + Q + R)^ + 

(p2 + m2)2 [^Q^ Q^B? (Q + i?)2 J + 4 p2 + „2 (^2 + ^2)2 2/^ ^4 

+^3 I ^ \ \ \ I 

Jpqr iP^ + ^2)2 g2 _|_ ^2 ,,2 ^ ^2 (p + q + r)2 + irfi ' 

Note that the second term is formally of order , but it vanishes at order e° due to the fact that Isnn = C'(e). 
The massive four-loop triangle sum-integral reads 

triangle(m ) - 2j.^^^^ g2 + ,^,2 (p + Q)2 + „j2 ^2 ^ ^2 (p + ^)2 + „j2 ^2 + ^2 (p + + ^2 ' (^-10) 

When all four momenta are hard, the leading contribution is given by setting m = 0, i. e. 

7-(hhhh) . 2n ^ v/ rr m 

When one of the momenta is hard and three are soft, we find 

t(1'ss«) frr)2l = ST^ / ^^——^ ^- \/ J_ + ... 

triangleV ) j^^^ p2 + ^2 ^2 ^ ^2 ^2 + ^2 (p + q + i.)2 + to2 ^Tg 5^ ^ 

This contribution is of order . When one momentum is soft and three momenta are hard, the contribution is 

1111 



2 



'pQR P2 + ^2 g2 + ^2 (p + Q)2 + ^2 ^2 + ^2 (p + p)2 + ^: 

This contribution is of order g"^ . When all four loop momenta are soft, the contribution is given by the massive 
three-dimensional triangle diagram /triangle: 



'"•^"g'*^^ ' ipg^ g2 + ^2 (p + q)2 + ^2 J.2 _^ ^2 (p + r)2 -H m2 s2 + ^2 (p -H s)2 -h m2 ■ ^' 

This contribution is of order g^ . Finally, we notice that the contribution when two momenta are soft and two momenta 
are hard, is of higher order in the coupling g. 



111111 
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Appendix D: EXPLICIT CALCULATIONS 

In this appendix, we illustrate the use of the calculational techniques developed by Arnold and Zhai in Ref. [l^ 
to evaluate complicated multiloop diagrams. The strategy is to rewrite the original sum-integral into two sets of 
terms. The first type is ultraviolet divergent, but is sufficiently simple to be evaluated analytically using dimensional 
regularization. The second type is finite both in the ultraviolet and the infrared, but is normally so complicated that 
it must be evaluated numerically. In order to isolate the divergences in terms that are tractable, typically one or more 
subtractions are required. 

We need to calculate the following three-loop diagram 

/-^^^{[n(P)P-(4^n(P)}, (D.i) 

where the self-energy n(P) is defined by 



rqQHP + Qr- 

The first term in Eq. l|D.ip arises from the m/T-expansion of the triangle sum-integral in four dimensions, while the 
second term arises from the term T/iXsun which is a part of the counterterm T2\,^\g^ I ■ 
At zero temperature, the self-energy is denoted by n"(P) and reads 

^ ' " (47r)2 [ ) r(2 - 2e) ' ^^'"^^ 

In order to isolate the UV divergences and simplify the calculations, we write the self-energy as 



n(P) = + n^iP) + n^(P) , (d.4) 



where n"(P) is the finite part of U°{P), i.e. we have subtracted the divergent piece in Eq. I|D.3P from 11° (P): 



and 11"^ (P) is the finite-temperature piece of n(P). In three dimension, n^(P) reads 







= ^/ (-"-f - I) . (D^6) 

where f = 2nTr. In the following we need the UV limit of 11"^ (P). This happens to be given by the UV limit of the 
full self-energy (|D.2p and is given by [l^ 



Using the decomposition (|D.4p . the integral in Eq. (|D.ip can be written as 



(4 



^ + j^mP)? + -^n°(p)n^(p) + ^[n^(P)]^ . (d.s) 
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We now consider the different contributions to /. The first term in Eq. I|D.8P is a simple one-loop sum-integral and 
reads: 



^ (47r)4e2 % p2 

2£ 7^2 



\AttTJ \2{AiiY 



^ + 2(1 



C'(-l)\ 1 TT 



C(-i) 



- + -T +4- 
e 4 



C(-i) C(-i) 



0(6) 



(D.9) 



The second term in Eq. IjP.sp contains no logarithmic UV divergences and so it is finite in dimensional regularization: 
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, C'(-i) 
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2 + log 



47rT 
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47rT 



0(6) 



(D.IO) 



The third term requires a little more thought. Since the UV behavior of n'^(P) is the integrand X\°^[P)\l^ [P) / 

is logarithmically divergent in the ultraviolet. In order to isolate this divergence, we add and subtract I\^y{P) from 
n2(P)n^(P)/P^. Thus the third sum- integral in Eq. p.sp becomes 

^3 = 2^ijn°(P)n^(P) 

= 2^ -^n°(p)[n^(p) - ul^iP)] + 2^ -LnO(p)n?;v(P) + ^n2(po - o,p)n^bo = o,p) ,(d.ii) 

where we have isolated the contribution from the po = term since the contribution to h from this term is infrared 
divergent. In order to calculate the first term in Eq. IjD.lip . we need Ii'^^{P) in coordinate space. It is given by the 
small-r behavior of 11^ (P) and reads 



nSv(^) 
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(47r) 
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This yields 



^3 = 2^^'-ijn2(P)[n^(P)~n5v(P)] 



IP 



d\- (^cothr---^ J2 



log 



2 1 ? 

Po+P 



The integral over 3-momentum can be done analytically. We write it as 



(2^)3p2+p2 



2 log 



(4^r)^ 
^°s;;2 



(D.13) 



(D.14) 



4ttT ^Pt,+P\ 

where the first two terms in the parentheses are independent of p, making this part of the integral a simple Fourier 
transform: 



d'^p 



(27r)3pg+p2 



2 + 2 log 



47rr 



Averaging over angles, the last term can be rewritten as 

d^p e*- 



{2T:fpl+p- 



log 2 I 2 
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log 2 , 2 
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(D.16) 
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The integrand has a branch cut starting aX p = i\po\ running to p ~ ioo, and a pole in p = «|po|- The contour can be 
deformed to wrap around the pole and the branch cut, and taking care to include contributions from both, one arrives 
at the result 

(Pp e'P" . f47rT)2 e-IPol"^ 



47rr 



, 2r 

log 7^ + 7B + e 
bo I 



2|po|r 



Ei(-2bok) 



where po = po/'^irT — n and the exponential-integral function Ei{z) is defined as 
Thus Eq. I|D.13|) can be rewritten as 



(D.17) 
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(D.19) 



The first three terms in the last parentheses are independent of r and po, and for these terms, the integral over r and 
the sum over Matsubara modes can be evaluated analytically. In particular, we are able to find the coefficient of log /i. 
This is fortunate, because it allows us to check the concistency of our final result for the free energy. Let 
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Integrating over angles and summing over Matsubara frequencies yields 
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The integral above is finite, but the individual terms are divergent for small r. We therefore regulate them by 
multiplying by an extra factor (2f)" and taking the limit a ^ in the end. The basic integrals we need are 

dtp 



e' - 
dtf" 



(e*-l)2 



r(x + l)C(a; + l), 

r(x + i) [c(x)-C(x + i)]. 



This yields 



AT' 
(4^ 



2r(a)[C(a - 1) - C(«)] + r(«)C(«) - 2r(a - 1)C(« - 1) - -r(a + l)C(a + 1) 

6 



The limit a — > is regular, and we obtain 



c = - 



2T2 
3(47r) 



1 + 7i; - 3 log(27r) + 2 
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c(-i)y 



(D.22) 
(D.23) 

(D.24) 

(D.25) 
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The remaining integral over the coordinate r as well as the Matsubara sum in Eq. (|D.19p must be done numerically. 
Eq. IjD.lQP can then be written as 



3 3(47r)4 
where the numerical constant C is 



(2 + 7£ + 21og^) h +7£;-31og(27r) + 2 
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C(-i) 



C 



c 



_3_ 

47r 



coth r 



1 f 

f ~ 3 



2f 



V e-2|po|riog _^_Ei(_2|pok) = 0.0034814, 
^ ^ IpoI ' 



The subtraction term in Eq. (|D.13p can be calculated with dimensional regularization and reads 
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The last term in Eq. (|D.lip is 

/g^ = 2T j 4n°(Po = o,p)n^(po = o,p) 

= 2T / 4n°(Po = o,p) [n(po - o,p) - n"(po = o,p)] , 

The second term vanishes in dimensional regularization since there is no mass scale in the integral, i.e. 

2T j ln°(po = o,p)nO(po = o,p) = 0. 



(D.26) 
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In order to evaluate the first term in Eq. l|D.29p . we must calculate n(po = 0,p). Using Feynman parameters, we 
obtain 
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(D.31) 



Inserting the expression for U^{po — 0,p) and Ii{po = 0,p), we obtain 
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where the prime indicates that we have omitted the po — mode from the sum. Expanding Eq. (|D.33P in powers of 
e, we obtain 
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The last term in Eq. p.8[l is 



h = 



5f ^[n^(^)]= 



(D.34) 



(D.35) 



Since the UV-behavior of U^{P) is l/P^, the sum- integral in Eq. I|D.35P is UV finite. However, I1t{P) has a loga- 
rithmic infrared divergence for the po = mode. This implies that the sum-integral h has linear and logarithmic IR 
divergences. The linear divergence is set to zero in dimensional regularization while the logarithmic is not. In order 
to isolate these divergences, we rewrite the sum-integral as 



(D.36) 



where the prime indicates that we have omitted the po = mode from the sum. The primed sum-integral in Eq. (|D.36P 
is finite both in the ultraviolet and in the infrared. Using the 3-dimensional representation of the IIt{P), Eq. (jD.6p . 
the first term in Eq. (jD.36P can be written as 
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The integral over three-momentum p corresponds to a Fourier transform of a massive propagator and so gives rise to 
a Yukawa potential. The sum over nonzero Matsubara frequencies can also be done analytically and we obtain 
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Averaging over angles, one finds 



/4 = jm / cothf - - cothf' - - 



(47r)4 Jo ''''' 
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log (e2(*^+^') - - log (^e^^+f^'+lf^-f^'l _ 1^ + |r _ f'l _ f _ f'j (D.39) 



The remaining integrals over r and r' must be done numerically and we obtain 

It - [0.0587392] . (D.40) 

The second term in Eq. l|D.36|) is rewritten as 

i"^ ^ T [ ^[n^(p„ = o,pr 

J p P 

^ T I ^ {[n^(po = o,p) - nfj,(p)]' + 2n^(po = o,p)nS,(p) - [iv(^{p)f} , (d.4i) 
J p p 

where Ili^{p) is given by the qo = term in Eq. l|D.3ip : 



n?p,(p) = T /^T-^^ : 



e-'-M'A 4^0? r(l/2 + e)r(l/2-e) 



= n^j (4.)3/^-^ r(i - .) — ■ 

The first integral in Eq. (|D.4ip is now well behaved in both the ultraviolet and the infrared. It can be evaluated 
numerically using the representation of n-^(po = 0,p) in three dimensions. The subtracted terms are infrared divergent 
and are calculated with dimensional regularization. The first integral can be calculated directly in three dimensions. 
In this case, Hf^ {p) reduces to 

^mip) = ^- (D.43) 
Using the three-dimensional representation (|D.6P for (P) with po = and Eq. I|D.43P , we get 
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The averages over the angles between p and r, and between p and r' can be done analytically and we obtain 



7-bl 



1 



1 



(4^)^ Jo Jo 
sin pr 



dr / dr' 



sm pr sm pr 
pr pr' 



coth r if coth r' 



dr- 



4(47r)p Jo pr 



coth r — — 
r 



1 



64p2 



(D.45) 



31 



The integrals over r, r', and p must be done numerically. The result is 

rp2 



li = ^[9.5763]. (D.46) 
The first subtraction term in Eq. (|D.4ip is 

J p P 

= 2T f ^ [Uipo = 0,p)- n"(po = 0,P)] Tlfj,{p) 

J p P 

^ 2T f ln(po = 0,p)nf^ip) , {D.47) 

J p p 

where we have used that the second term vanishes in dimensional regularization: This term is logarithmically divergent 
both in the infrared and in the ultraviolet. If we use the same scale for the regularization of ultraviolet and and infrared 

Inserting the expressions for n°(po = 0,p) and nfj^(p) into Eq. l(DA7\i . we obtain 
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_ / e^^^i^ \ 4^ r(i + e)r(l - e)r(-2e)r(l + 3e)r2(l + 2e) 

~(4^1v4^2^j 7 r(|-e)r(l-e)r(2 + 4e) 

The prime on the sum in the second line indicates that we have excluded the zero mode, qo = 0, from the sum. This 
mode gives rise to an integral that is linearly divergent in the infrared. Since there is no mass scale in this integral, it 
vanishes. Note also that the integral over p is logarithmically divergent in the infrared and this divergence is not set 
to zero in dimensional regularization jl^l- Expanding Eq. I|D.48P in powers of e, we obtain 
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+367_E log -/^ - 127b - 3671 



Finally, we consider the last subtraction term in Eq. (|D.4ip . Since 11^ (po = 0,p) goes like 1/p for small p, the 
integrand has a linear infrared divergence. This divergence is set to zero in dimensional regularization. In fact, since 
there is no mass scale in the integral, it vanishes: 

T I 4[n?a(p)]^ = 0. (D.50) 

Jp P 

Adding Eqs. ijPlOl) . (|D361) . (|D381) . |(D34| . (|R401) . (|D^ . and (|R49l) . we can write / in the following form: 
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